CePtaSi: an unconventional superconductor without inversion center 
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Most superconducting materials have an inversion center in their crystal lattices. One of few 
exceptions is the recently discovered heavy-fermion superconductor CePtsSi [E. Bauer et al, Phys. 
Rev. Lett. 92, 027003 (2004)]. In this paper, we analyze the implications of the lack of inversion 
symmetry for the superconducting pairing. We show that the order parameter is an odd function of 
momentum, and that there always are lines of zeros in the excitation energy gap for one-component 
order parameters, which seems to agree with the experimental data. The superconducting phase can 
be non-uniform, even without external magnetic field, due to the presence of unusual gradient terms 
in the Ginzburg-Landau free energy. Also, we performed ab initio electronic structure calculations 
for CePtsSi, which showed that the spin-orbit coupling in this material is strong, and the degeneracy 
of the bands is lifted everywhere except along some high symmetry lines in the Brillouin zone. 



PACS numbers: 74.20.Rp, 74.70.Tx, 71.20.Be 



I. INTRODUCTION 

In the past two decades, a number of novel supercon- 
ducting materials have been discovered where order pa- 
rameter symmetries are different from an s-wave spin sin- 
glet, predicted by the Bardeen-Cooper-Schrieffer (BCS) 
theory of electron-phonon mediated pairing. From the 
initial discoveries of unconventional superconductivity in 
heavy-fermion compounds, the list of examples has now 
grown to include the high- Tc cuprate superconductors, 
ruthenates, ferromagnetic superconductors, and possibly 
organic materials. In most of these materials, there are 
strong indications that the pairing is caused by the elec- 
tron correlations, in contrast to conventional supercon- 
ductors such as Pb, Nb, etc. Non-phononic mechanisms 
of pairing are believed to favor a non-trivial spin struc- 
ture and orbital symmetry of the Cooper pairs. For ex- 
ample, the order parameter in the high- Tc superconduc- 
tors, where the pairing is thought to be caused by the 
anti- ferromagnetic correlations, has the d-wave symme- 
try with lines of zeroes at the Fermi surface. 

A powerful tool of studying unconventional supercon- 
ducting states is symmetry analysis, which works even 
if the pairing mechanism is not known. Mathemati- 
cally, superconductivity is a consequence of the break- 
ing of the gauge symmetry U{1) in the full symmetry 
group of the normal state: Q = Sm x U{1), where Sm 
is a magnetic space group, which includes usual space 
group operations (i.e. translations, rotations, inversion, 
etc) and time reversal operation K. In a non-magnetic 
crystal, Sm = S x JC, where S is the space group, and 
K. = {E, K}. In a magnetic crystal, K enters S]\i only 
in combination with other symmetry elements. The su- 
perconducting state is said to be "unconventional" if, in 
addition to the gauge symmetry, some other symmetries 
from Q are broken. The group-theoretical analysis of 
unconventional superconducting states in non-magnetic 
crystals was developed in Refs. Recently jit was 

extended to include the ferromagnetic case 

In almost all previous studies, it was assumed that the 



crystal has an inversion center, which leads to degener- 
ate electron bands and makes it possible to classify the 
Cooper pair states according to their spin (or pseudospin 
if the spin-orbit coupling is taken into account). The 
even (singlet) and odd (triplet) components of the su- 
perconducting order parameter can then be studied sep- 
arately 0. Although this is the case in most supercon- 
ductors, there are some exceptions. Early discussion of 
the possible loss of inversion symmetry associated with a 
structural phase transition in VsSi, which is an A15 type 
superconductor, can be found in Rcf. Later, a C15 
superconductor HfV 2 was found to undergo a transition 
from a cubic to a non-centrosymmetric body-centered or- 
thorhombic structure f^. The possible existence of su- 
perconductivity was reported in ferroelectric perovskite 
compounds SrTiOa and BaPbOs-BaBiOs [13. More 
recently, there has been a renewed theoretical interest in 
this problem, mainly in the context of superconductivity 
in two-dimensional (2D) systems, such as Cu-0 layers in 
YBCO [ijl, or surface superconductivity on Tamm lev- 
els 0, Il3l |. According to Ref. 0, , in the absence of 
inversion symmetry the order parameter becomes a mix- 
ture of spin-singlet and spin-triplet components, which 
leads, for instance, to the Knight shift attaining a non- 
zero value at T = 0. Other peculiar properties of non- 
centrosymmetric superconductors include a jump of mag- 
netic induction at the surface and the possibility of 
non-uniform helical superconducting phases due to the 
presence of first-order gradient terms in the Ginzburg- 
Landau (GL) functional TB'l. 

This work is motivated by a recent experimental dis- 
covery of superconductivity with Tc ~ 0.75K in CePtsSi 
|l6l |. which is the first known heavy-fermion material 
without inversion center. (It should be mentioned that 
incommensurate density modulations might break inver- 
sion symmetry in a heavy-fermion compound UPta, as 
pointed out in Ref. 17].) Our goal is two-fold. First, we 
study the symmetry of electron bands, calculate the elec- 
tronic structure, and estimate the magnitude of the spin- 
orbit coupling in the normal state of CePtaSi. Second, we 
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use general symmetry arguments to analyze the gap sym- 
metry and spatial dependence of the order parameter in 
a three-dimensional non-centrosymmetric tetragonal su- 
perconductor, assuming a strong spin-orbit coupling and 
the clean limit. The paper is organized as follows. In 
Sec. m we study symmetry of the single electron states. 
In Sec. mil the results of the electronic band structure 
calculations are reported. Sec. IIVI focuses on the sym- 
metry of the superconducting order parameter, possible 
locations of the gap nodes, and also spatial structure of 
the superconducting phase, using the phenomenological 
Ginzburg-Landau approach. Sec. concludes with a 
discussion of our results. 



II. SYMMETRY OF ELECTRON BANDS 

The symmetry group of the normal paramagnetic 
state: Q = S x IC x U{1), where S is the space group 
of the crystal, and C/(l) is the gauge group. In the case 
of CePtsSi, the space group is P4mm (No. 99), which 
is generated by (i) lattice translations by the primitive 
vectors a = a(l,0, 0), b — a(0, 1,0), c — c(0,0, 1) of a 
tetragonal lattice, and (ii) the generators of the point 
group G = C^v- the rotations C^z about the z axis by 
an angle 7r/2 and the reflections cFx in the vertical plane 
(100). Spatial inversion / is not an element of the sym- 
metry group. 

In the absence of inversion symmetry, spin-orbit (SO) 
coupling plays a crucial role due to its band-splitting ef- 
fect. At non-zero SO coupling, the single-particle wave- 
functions are linear combinations of the eigenstates of 
the spin operator s^: (r|7/>) = u(r)x'\ +v{r)xi- Since the 
normal state Hamiltonian Hq is invariant with respect to 
the crystal lattice translations, the eigenfunctions are the 
Bloch spinors 4'k,n{f) belonging to the wave vectors k in 
the Brillouin zone, which can be written in the form 

(r|fc,n) = 'Ufe,„(r)xT +Vk,ni'^)xi- (1) 

The corresponding eigenvalues e„(fc) describe the band 
dispersion of free electrons. 

In the presence of both the time-reversal and inversion 
symmetries the bands are two-fold degenerate at each k. 
Indeed, the states \k,n) and KI\k,n) correspond to the 
same fc, have the same energy, and are orthogonal (in 
addition, these two states are degenerate with another 
pair of orthogonal states, K\k, n) and /|fc, n), which cor- 
respond to — fc). In this case, n — {v, ±), where ± labels 
the linearly independent Bloch states at a given band 
index v. There is a freedom in choosing the basis func- 
tions \k,iy,+) and — ). The most frequently used 
convention is that they should transform under rotations 
R similar to the spin eigenstates \k, v, f) and |fc, J,), i.e. 

R\k,v,a) ^Uc,l3{R)\Rk,v,f3), (2) 

where a,/3 = ± and U{R) is the spin-rotation matrix: 
for a rotation by an angle 9 around some axis n: U (i?) = 



exp[— i(6'/2)(cr • n)] (cr are Pauli matrices) . The practical 
recipe for constructing the basis of the Bloch states is 
as follows: first choose a state ly, +) at each k in the 
irreducible part of the first Brillouin zone, then act on it 
by KI to obtain an orthogonal state \k,iy,~) at the same 
fe, and finally act on |fc, i^, ±) by the elements of the point 
group and use the prescription ^ to obtain the pairs 
of the basis functions belonging to the star of k. The 
single-electron states constructed in this way are referred 
to as the pseudospin states |^]. It is the presence of the 
inversion center that makes the bands double degenerate 
in a non-magnetic crystal: e,y^+{k) = e,y,+ (— fc) = e^,-{k) 
at all k. 

In contrast, if the crystal lacks the inversion symmetry, 
then the degeneracy of the single-electron bands en{k) is 
lifted everywhere, except from some points or lines of 
high symmetry. In particular, the bands always touch at 
fe = (the r point), because of the Kramers theorem: 
time-reversal symmetry means that \k,n) and K\k^n) ~ 
— k,n) have the same energy. In the case of CePtaSi, 
one can show that the spin-split bands touch along the 
rZ line (i.e. along the [001] direction), and there are 
no other symmetry-imposed degeneracies in the bands 
crossing the Fermi level (see Sec. IIIII below). In the limit 
of zero SO coupling (which is not applicable to CePtaSi), 
but still without an inversion center, the symmetry of the 
system contains SU (2) spin rotations, in addition to the 
space group, time reversal, and the gauge group. Then, 
the bands at each k are double degenerate. 

Let us now show that the transformation properties 
of the single-electron wave functions in the absence of 
inversion center are different from those of the spin or 
pseudospin eigenstates, see Eq. 0. Mathematically, 
the wave functions in each band transform according 
to irreducible co-representations of the Type II mag- 
netic space group Sm — S x K. 0| (one should use 
co-representations because the time-reversal operation 
K is anti- unitary). In addition, the co-representations 
must be double- valued because the states |fe) are spin-1/2 
spinors, so that any rotation by 27r changes their sign: 
Cl,\k) ^ -|fe), al\k) = -|fe), and also K'^\k) = -\k) 
(the band index n is omitted). The double- valuedness 
can be dealt with in the standard fashion by using the 
"double-group" trick : one introduces a fictitious new 
symmetry element which commutes with all other 
elements and satisfies the conditions Cl^ = a'^ = E, 
E^ = E, and also K"^ = E. 

The co-representations of Sm can be derived from the 
usual representations of the unitary component, which in 
our case coincides with the space group S. For each fe in 
the Brillouin zone, the basis of an irreducible represen- 
tation of S is formed by the Bloch states corresponding 
to the star of fe. The state K\k) belongs to the wave 
vector — fe, but the irreducible representations of S cor- 
responding to fe and — fe are inequivalent (because of the 
absence of inversion symmetry) and must therefore be 
regarded as belonging to a single "physically irreducible" 
representation of twice the dimension 19] . In terms of co- 
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representations this means that the Bloch states |fe) and 
K\k) belong to the same two-dimensional irreducible co- 
representation of S* X /C 18]. Thus, the appropriate basis 
of the Bloch states is formed by the wave functions \Gk) 
corresponding to the star of fc, and also by their time- 
reversed counterparts K\Gk), which can be combined in 
a set of bispinors l^'fc) = {\k), K\k))'^ . All the states from 
this set have the same energy: e(fc) = e{Gk) — e{—k). 

Since the function G\k) belongs to the wave vector 
Gk, one can write G\k) = e^'^'^'-'^^Gk) . The undeter- 
mined phase factors come from the freedom in choos- 
ing the phases of the Bloch states and realize a rep- 
resentation of the point group in this basis (it is as- 
sumed that the Bloch states are single-valued and con- 
tinuous functions throughout the Brillouin zone). Us- 
ing the commutation of G and K we have GK\k) — 
KG\k) = K ey.Y>[iVk{G)]\Gk) = ey.^[-iipk{G)]K\Gk) . 
The co-representation matrices in the basis of bispinor 
wave functions can be obtained from the relations 

/ „i(pk{G) n \ 

Gi*.) = (' j (3) 
m>.) = ( V ) i^fc) (4) 

^I*fc)=(_"l (5) 

[we have taken into account that K{K\k)) — — |A;)]. Note 
that the multiplication rules for the co-representation 
matrices are different from usual unitary group repre- 
sentations. For example, D(GiG2) = D{Gi)D{G2) and 
D{GK) ^ D{G)D{K), but D{KG) = D{K)D*{G) and 
dIk"^) = D{K)D*{K) tl8|. 

The general symmetry arguments given above can be 
illustrated using an exactly-solvable three-dimensional 
generalization of the Rashba model, which was originally 
proposed to describe the effects of symmetry lowering 
near the surface of a semiconductor '20l| and recently ap- 
plied in Refs. [H 113 to quasi-2D non-centrosymmetric 
superconductors. Consider a single band e{){k) in a crys- 
tal described by the point group C^y. At zero SO cou- 
pling, the band is two-fold degenerate due to spin. The 
absence of reflection symmetry in the xy plane implies 
the existence of an internal electric field in the crystal, 
whose average over a unit cell is non-zero. In the Rashba 
approximation, this non-uniform field is replaced by its 
average, introducing a constant vector n \\ z. When a 
non-zero SO coupling is switched on, it can described by 
an additional term in the Hamiltonian: 

Hso ^ a'^n ■ {(T„„, y. k) a\^ak,a', (6) 
fe 

where cr, cr' —^,[ is the z-axis spin projection, and the 
states |fe,cr) are the Bloch spinors at zero SO coupling. 
Diagonalization of the full single-electron Hamiltonian, 
H — Hq + Hso, gives two bands 

ei(fc) = eo(fc) +a|fci|, e2{k) ^ eo{k) ~ a\k±\ (7) 



{\k±\ = + fc^), which satisfy the condition ei.2(fc) = 

^i,2{—k) and additional symmetries from the point 
group. Also, the bands touch along the line k \\ z. It 
is easy to see that one cannot use pseudospin to label 
these bands, because the time reversal K transforms the 
bands into themselves: n) ^ \ —k,n) [n — 1, 2), while 
the pseudospin states would be transformed into one an- 
other. 

To conclude this section, it should be mentioned that 
the superconductivity in CePtaSi seems to occur in the 
presence of anti-ferromagnetic (AFM) order al- 
though no data have been reported on the structure of the 
magnetic phase or the magnitude of the staggered mo- 
ment. Although the symmetry analysis above was done 
assuming a paramagnetic normal state, our results can 
be easily generalized for an AFM case. For a staggered 
magnetization directed along the z axis, the only change 
one has to make in the symmetry group is to replace K 
with KTa, which combines the time reversal operation 
with a lattice translation. Then, Eq. Q is replaced by 

/ Q ika \ 

KT^\'^u)=^_^-,Ua j l^fe), (8) 

so that {KTa.f\k) = -e-2«*=°|fc). 

From the expressions Q-®, one obtains the trans- 
formation rules for the creation operators of electrons in 
the Bloch states |fc). Although there is some freedom in 
choosing the phase factors, we will see in Sec. |^ that 
the physically relevant properties are insensitive to the 
choice of (pk{G). 

III. ELECTRONIC STRUCTURE 

CePtaSi crystallizes in the same tetragonal structure as 
CePtsB, with space group P4mm (No. 99). The lattice 
parameters are a = 4.072 A and c = 5.442 A. Ce is at 
the 1(b) site (1/2, 1/2, 0.1468), Pt at the 2(c) site (1/2, 
0, 0.6504) and at the 1(a) site (0, 0, 0). The z coordinate 
of the latter site was chosen to be zero. Si is at the 1(a) 
site (0, 0, 0.4118). These structural parameters all derive 
from single-crystal x-ray data |T^ . and can be assumed 
to be sufficiently accurate for our purposes. 

We calculated the electronic structure of CePtsSi 
with the full-potential (FP) linear augmented-plane wave 
(LAPW) method l2lL 122. l23| . which is based on density 
functional theory |2J|. For the exchange and correla- 
tion potential we used the local density approximation 
[25. 26] (LDA). We performed a non-magnetic calcula- 
tion, neglecting the AFM order observed experimentally 
below 2.2 K [l^. The muffin-tin radii of the atoms were 
chosen as 2.11 oq. A typical number of plane waves in 
our basis set was 580. The electronic ground state was 
calculated self-consistently on a grid of 4212 k points in 
the entire Brillouin zone. 

For the electronic structure of alloys containing heavy 
elements, such as Ce and Pt, the SO coupling can in 
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general not be neglected. In particular, as shown in Sec. 
IIVI for the analysis of the symmetry properties of the 
superconducting order parameter it is important to ob- 
tain an estimate of the SO splitting of the bands near the 
Fermi energy. The SO coupling has been included in our 
calculations using the "second variational treatment" , as 
discussed by MacDonald et al [l^. In this approach, first 
the eigenstates are calculated in the absence of the SO 
interaction. Then, the SO interaction is included in a 
perturbative way, where the eigenstates up to a certain 
cut-off energy, calculated without the SO interaction, are 
used as the basis states. In our calculations this cut-off 
energy was approximately 22 eV above the Fermi energy. 

Although the CePtsSi structure doesn't have inver- 
sion symmetry, the band energies still satisfy the rela- 
tion e„(k) — £„(— k), due to the time reversal symmetry 
of the single-electron Hamiltonian (see Sec. EJ. Figure 
n shows the band structure of CePtaSi without SO cou- 
pling calculated along some high symmetry lines. The 
bands were plotted according to increasing energy. A 
complete analysis of band crossings based on the charac- 
ter of eigenfunctions was not carried out. However, the 
bands labeled (3 and 7 do cross between the symmetry 
points X and M, and so do the (3 and 7 bands between 
R and A. The 7 band and the first dotted band above 
the Fermi energy have the same symmetry between M 
and r, and hence are unlikely to cross there. The labels 
of the bands were chosen according to the band index, 
not the band character, simply to relate various parts of 
the Fermi surface to the bands crossing the Fermi level. 
In the electronic density of states (not shown) there is a 
peak at 0.4 eV above the Fermi energy, which is due to 
the unfilled Ce-4f electrons. We found that the electrons 
at the Fermi energy are predominantly of Ce-4f charac- 
ter. 

Figure [3 shows the band structure of CePtaSi with SO 
coupling (22|]. As in Fig. ^we connected bands with the 
same band index, ignoring band crossings. In Fig. [3 the 
bands near the Fermi energy tp are split by an amount 
of at most 50 — 200 meV. This splitting vanishes along 
the lines T - Z and M - A. 

FigureOshows the cross sections of the Fermi surface of 
CePtsSi in the presence of SO couphng. The splitting of 
the bands at the Fermi energy due to the SO coupling is 
prominent in this figure. The sheets of the Fermi surface 
labeled a and /? are hole-like. The 7 sheets are electron- 
like. The a sheet consists of a feature around the Z point. 
The [3 band gives rise to a larger feature around the Z 
point and, in addition, to a feature around the X point 
[and the symmetry related point Y = (0,0.5,0)]. The 7 
bands give rise to a feature around the line M — A, which 
is almost dispersionless in kz , as well as to a small feature 
around the F point. We further noted that the a bands 
contribute 1.9% and 3.5% respectively to the density of 
states at the Fermi energy. Similar contributions coming 
from the SO split (3 and 7 bands are 25% and 45% and 
15% and 9.0%, respectively. 



IV. SUPERCONDUCTING ORDER 
PARAMETER 

A. Symmetry analysis 

The single-electron states |fc, n) can be used as a basis 
for constructing the Hamiltonian which takes into ac- 
count the Cooper pairing between electrons. We have 



Hsc, where 



(9) 



is the free-electron part, with the chemical potential ab- 
sorbed into the band energies. As follows from the results 
of Sees. IhI and HTTl the electronic bands e„(fc) are non- 
degenerate, except along the line k \\ z, and invariant 
under all operations from the point group C^y and also 
under inversion: en{k) = e„(— fe). The Fermi surface of 
CePtsSi consists of several sheets (see Sec. Illljl . all of 
which can in principle participate in the formation of the 
superconducting order. 

Assuming a BCS-type mechanism of pairing, the inter- 
action between the band electrons in the Cooper channel 
can be written in the following form: 



(10) 



where 



^ic' = jEE^"'nfc,fc')4.„cl.,„^-^',nC.',„ (11) 



n k,k' 



^ic^ - \ E E^«™(^'fc')4,„cl,_„c_fe,,„c,,™ (12) 



n^m h.h' 



Hif = \ E E^»-(^''^')4,„cL,,,„c_fc,,„.Cfe,.„. (13) 

n^m h.h' 

The potentials V are non-zero only inside the energy 
shells of width lJc (the cutoff energy) near the Fermi sur- 
face. 

Treating the Cooper interaction between the electrons 
with opposite momenta in the mean-field approximation, 
we obtain 



Hsc — 



^EE[^«"('^)4,,/-fc,,„+h.. 



(14) 



h nni 



Here the diagonal matrix elements A„„(fe) represent 
the Cooper pairs composed of quasiparticles from the 
same sheet of the Fermi surface, and the off-diagonal 
matrix elements A„m(fc) with n ^ m represent the 
pairs of quasiparticles from different sheets. From the 
anti-commutation of the fermionic operators, it follows 
that A„„(A;) are odd functions of k, but A„m(fc) = 
— Am„(— fc) with n ^ m do not have a definite parity. 
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A considerable simplification occurs if to assume that 
the superconducting gaps are much smaller than the in- 
terband energies. The band structure calculations of Sec. 
mil show that typically the SO band splitting Ego is of 
the order of 50 — 200 meV (between the bands derived 
from the degenerate spin-up and spin-down bands at zero 
SO coupling), which exceeds the superconducting gap 
by orders of magnitude. In this situation, the forma- 
tion of interband pairs is energetically unfavorable, and 
the amplitudes A^m with n ^ m vanish. The origin 
of the suppression of these types of pairing is similar to 
the well-known paramagnetic limit of singlet supercon- 
ductivity |23- the interband splitting Eso cuts off the 
logarithmic singularity in the Cooper channel, thus re- 
ducing the critical temperature. Although the condition 
Ego 3> Tc is violated very close to the poles of the Fermi 
surface where the spin-split bands touch, the off-diagonal 
Cooper pairing in the vicinity of these points is still sup- 
pressed due to the phase space limitations. We also ne- 
glect the possibility of the Cooper pairs having a non-zero 
momentum, i.e. (c^^g n'^-fc m) ^ ^ [Larkin-Ovchinnikov- 
Fulde-Ferrell (LOFF) phase] Although the critical 

temperature of the resulting non-uniform superconduct- 
ing state can be higher than that of the uniform state, 
this would not be sufficient to overcome a large depairing 
effect of the SO band splitting. 

Thus, the interband pairing H13|l can be neglected, and 
A„„i(fc) = An{k)5mn [sfl • This is reminiscent of the sit- 
uation in ferromagnetic superconductors, where only the 
same-spin components of A survive the large exchange 
band splitting Q. It follows from Eqs. ®-© and anti- 
linearity of K that 

Eici . )E-^ = cl . 

V k.n —k,7i/ k,n — fe.n 

where A is an arbitrary constant. We have the following 
transformation rules for the order parameter under the 
elements of Q: 

G : A„(fc) ^ A„(G-ifc), 

K : A„(fc)^A;(fc). ^^^^ 

Thus, the order parameter components transform like 
scalar functions. There is no need for double groups, 
since E is equivalent to E when acting on A. In the 
case of an AFM normal state, K should be replaced with 

The superconducting order parameter on the nth sheet 
of the Fermi surface transforms according to one of the 
irreducible representations F of the normal state point 
group Civ It can be represented in the form 

A„,r(fc) ^ rin,i4'r,n,i{k), (16) 

i=l 



where dr is the dimension of F, (j)T,n,i{k) are the basis 
functions (which are different on different sheets of the 
Fermi surface in general), and rjn^i are the order param- 
eter components that enter, e^., the GL free energy and 
can depend on coordinates Despite the absence of 
inversion center in the crystal, the order parameters A„ 
have a definite parity, namely they are all odd with re- 
spect to fc — > — fe. The odd irreducible representations of 
^-"41; are listed in Table IJ Since dp = 1 or 2, the order 
parameter in each band can have one or two components. 

(2) 

If we neglect the interband pairing described by Hsc , 
see Eq. (|12f) . then the order parameters A„ are com- 
pletely decoupled, in particular, they all have different 
critical temperatures Tc^„. However, there is no reason 
to expect these interband terms to be small. If they are 
taken into account, then all A„(A;) are non-zero, so the 
superconductivity will be induced simultaneously on all 
sheets of the Fermi surface. The simplest way to see how 
this works is to use the GL free energy functional, which 
contains all possible uniform and gradient terms invariant 
with respect to Q. For a one-dimensional representation 
F (the generalization for two-dimensional representations 
is straightforward), we obtain the following expression for 
the uniform terms in the free energy density: 

Funiform = ^ A„miT)ri^r]„i + F4, (17) 



where stands for the terms of fourth order (and 
higher), and Aij is a real symmetric matrix. The off- 
diagonal elements of A correspond to the interband pair- 
ing. The critical temperature Tc is defined as the maxi- 
mum temperature at which A ceases to be positive defi- 
nite. Below Tc all rjn are non-zero and proportional to a 
single complex number rj, such that rjn = enrj, where £„ 
are constants that can be found by minimizing Funiform- 
Therefore, the components of a one-dimensional order 
parameter corresponding to the representation F are 
given by 

A„(fc) (fe). (18) 

Although the basis functions (f) are different in general, 
they all have the same symmetry. 

Our phenomenological theory cannot determine which 
pairing channel corresponds to the highest critical tem- 
perature. In a recent work, Frigeri et al 'si'l proposed 
a microscopic model for superconductivity in a non- 
centrosymmetric crystal, treating the SO coupling as a 
perturbation with a Rashba-type Hamiltonian. Assum- 
ing a strong interband pairing interaction which induces 
order parameters of the same amplitude on both sheets of 
the Fermi surface [see Eq. (U))], they predicted a certain 
gap symmetry for CePtaSi, which seems to correspond 
to the two-dimensional E representation in our classifi- 
cation. 
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TABLE I: The character table and the examples of the odd 
basis functions for the irreducible representations of C41, . 



r 


E 


C4. 




'^r(fc) 


Ai 


1 


1 


1 




A2 


1 


1 


-1 


{kx ky^kxkyhz 


Bi 


1 


-1 


1 




B2 


1 


-1 


-1 


kx ky kz 


E 


2 








kx J ky 



B. Gap zeros 

The symmetry considerations can help find the zeros in 
the energy spectrum of Bogoliubov quasiparticles, which 
are responsible for pecuharities in the low-temperature 
behavior of unconventional superconductors . The gap 
structure of the order parameter belonging to the two- 
dimensional representation E of C41, depends on the su- 
perconducting phase, i.e. on the values of the compo- 
nents r]i and 772, which in turn are determined by mini- 
mizing the free energy of the superconductor. In contrast, 
the gap nodes for the one-dimensional order parameters 
are required by symmetry. Although the momentum de- 
pendence of the order parameter is different on different 
sheets of the Fermi surface, see Eq. (jlSII . the locations of 
symmetry-imposed gap zeros are the same. 

One can easily show that the order parameter corre- 
sponding to Ai vanishes on the plane kz = 0, so that the 
energy gap has lines of nodes at the equators of all sheets 
of the Fermi surface. Indeed, 



(19) 



therefore, (pAiikx, ky,0) = 0. In a similar fashion, one 
can prove the existence of lines of zeros a.t kz = for all 
other one-dimensional representations. Also, the basis 
functions of A2 and B2 have lines of zeros a.t k^ — 
or ky = 0, while the basis functions of A2 and Bi have 
lines of zeros at kx = iky. The examples of the basis 
functions that have only zeros imposed by symmetry are 
given in Tabled 

From the results of Sec. IIIII it follows that some of 
the sheets of the Fermi surface cross the boundaries of 
the Brillouin zone. As seen from Eq. H19|) . the order pa- 
rameters corresponding to all one-dimensional represen- 
tations vanish at kz = ±7r/c, i.e. at the top and bottom 
surfaces of the Brillouin zone, because {kx,ky,TT/c) and 
{kx, ky, —tt/c) are equivalent points. In addition, for the 
same reason the basis functions of A2 and B2 have lines 
of zeros at kx = iir/a or ky = ±7r/a. 

The gap nodes for one-dimensional order parameters 
are present at the same locations on all sheets of the 
Fermi surface and can disappear only if the interband 
pairing interactions Hi^\ see Ec^. (|13|) . are taken into 
account. 



C. Helical superconducting states 

In addition to the uniform terms (|17() . the GL func- 
tional for the order parameter corresponding to a one- 
dimensional representation of C^y also contains gradient 
terms 

Fgrad,l ^^K^^{D±rin)*{Dj_r].^) 



+ KmiDzVn)*iDzr]m) 



(20) 



where D = V + i(27r/$o)A, <i>o = nhc/e is the flux 
quantum, A is the vector potential, and D± = {Dx, Dy). 
The coefficients and are real symmetric ma- 

trices. The terms (|2U|I . which are of the second order 
in D, are present in any multi-band tetragonal super- 
conductor with the order parameter corresponding to a 
one-dimensional representation of the symmetry group. 
However, in the absence of an inversion center, one can 
have additional terms in the GL functional, which sat- 
isfy all the necessary symmetry requirements but are of 
the first order in gradients In our case, they can be 
written in the form 

Fgrad,2 = LnniiVnDzVm " TJ^DzTIn), (21) 



where Lnm is a real anti-symmetric matrix, which is non- 
zero only if the interband pairing H12I) is present. The 
terms (I21II lead to the possibility that the superconduct- 
ing state which appears immediately below Tc can be 
non-uniform, even without external magnetic field. 

Consider for simplicity only two bands participating 
in superconductivity, i.e. n = 1,2. In this case, the 
matrix L in Eq. (|21() has only one non-zero element: 
Li2 = —L21 = A/2. The critical temperature for the 
superconducting state 



(22) 



is determined by the stability condition of the quadratic 
terms (both uniform and gradient) in the GL functional 
towards formation of a state with non-zero rin,o- From 
Eqs. H17I20I21|I . one has the following equation for Tc{q): 



det 



An + Kiiq^ A12 + Kuq^ + iXq 
A12 + Ki2q^ - iXq A22 + K22q^ 



= 0, (23) 



where Aii(T) = ai(T - Ti) and A22{T) = a2{T -T2), 
with Ti and T2 having the meaning of the critical temper- 
atures for the bands 1 and 2 respectively in the absence 
of any interband coupling. The phase transition tem- 
perature is obtained by the maximization of Tc{q) with 
respect to q. It is easy to show that the maximum crit- 
ical temperature corresponds to a state with g ^ (a 
"helical" state ^3) if the following condition is satisfied: 

A2 + 2A12K12 - An{T,,^)K22 - A22{Tcfi)Kn > 0, (24) 
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where Tc_o = Tc{q = 0). However, even if this condi- 
tion is violated and the phase transition occurs from the 
normal state to a uniform superconducting state, there 
remains a possibility that this uniform state becomes un- 
stable towards the formation of a helical state at a lower 
temperature. To find this instability, one has to include 
non-linear terms in the free energy. Because of a large 
number of the phenomenological parameters in the GL 
functional with the higher-order terms, the phase dia- 
gram of this system is quite rich |l5j . In particular, there 
exist various types of helical phases with g 7^ 0, separated 
from one another and from the uniform phase by addi- 
tional phase transitions below Tc. 

It should be emphasized that the origin of the heli- 
cal superconducting states is different from that of the 
LOFF non-uniform states . In terms of the GL func- 
tional, the LOFF state corresponds to the sign change 
of the second-order gradient term at some values of the 
parameters (e.g. of the external magnetic field), while 
our helical instability occurs because of the presence of 
the first-order gradient terms. 

V. CONCLUSION 

We have shown that the order parameter in a non- 
centrosymmetric superconductor with strong spin-orbit 
coupling has only intra-band components and is always 
odd with respect to k —k, which is a consequence of 
Pauli principle. This should be contrasted to the case 
of zero spin-orbit coupling, in which the bands are two- 
fold degenerate. In that limit, one cannot separate the 
odd and even components of A because of the lack of 
inversion symmetry, so the order parameter does not have 
a definite parity. 

The Fermi surface of CePtaSi consists of three pairs of 
sheets, a, (3, and 7, each split by the spin-orbit coupling. 
Our band structure calculations reveal that the states 
at the Fermi level are predominantly of Ce-4/ character. 
These states are affected strongly by spin-orbit coupling, 
which leads to the band splitting energy as high as 50 — 
200 meV. The splitting vanishes along the T - Z and A 
- M symmetry lines. By far the biggest contribution to 
the density of states at the Fermi level comes from the (3 
sheets. 

Although the large value of the spin-orbit band split- 
ting excludes the superconducting states that correspond 
to the pairing of electrons from different sheets of the 



Fermi surface, one can expect that the interband pair 
scattering will induce gaps of the same symmetry on all 
sheets of the Fermi surface. Possible gap structure of 
CePtsSi depends on the dimensionality of the order pa- 
rameter. If the order parameter corresponds to a one- 
dimensional representation of the group C4i,, then the 
gap has line nodes where the Fermi surface crosses the 
high-symmetry planes or the boundaries of the Brillouin 
zone: at kz — 0, ±7r/c for all ID representations; at 
kx — 0, ±7r/a and ky — 0, ivr/a for A2 and B2', at 
kx — iky for A2 and Bi. 

The presence of the gap nodes would manifest itself, 
e.g. in a power-law behavior of thermodynamic and 
transport characteristics at T ^ 0. Although the gap 
symmetries on all sheets of the Fermi surface should be 
the same, their magnitudes may be different. The exper- 
imental dataj_e.g. a reduced value of the specific heat 
jump at Tc 16], indicate that only some parts of the 
Fermi surface have non-zero superconducting gaps, while 
others remain normal. If this is indeed the case, then 
the specific heat would drop as G{T)/T cx const + aT at 
low temperatures (with the constant contribution com- 
ing from the normal sheets of the Fermi surface), which 
seems to agree with the experimental data of Ref. 16] in 
zero field. More detailed information about the pairing 
symmetry can be obtained only if the precise location of 
the line nodes is known. 

The absence of inversion symmetry can also have in- 
teresting consequences for the spatial structure of the 
superconducting phase. We showed that, in contrast 
to the centrosymmetric case, the Ginzburg-Landau free 
energy can now contain additional terms which are of 
the first order in gradients. Such terms can make the 
superconducting phase unstable towards the formation 
of a non-uniform (helical) state even at zero magnetic 
field. For the order parameters corresponding to the one- 
dimensional representations of C41, , this possibility exists 
only if the interband pairing is taken into account. 



Acknowledgements 

The authors thank B. Mitrovic for useful discussions, 
and also D. Agterberg, E. Bauer, and especially V. Mi- 
neev for stimulating correspondence. This work was sup- 
ported by the Natural Sciences and Engineering Research 
Council of Canada. 



[1] G. E. Volovik and L. P. Gor'kov, Sov. Phys. JETP 61, 

843 (1985) [Zh. Eksp. Teor. Fiz. 88, 1412 (1985)]. 
[2] K. Ueda and T. M. Rice, Phys. Rev. B 31, 7114 (1985). 
[3] K. V. Samokhin and M. B. Walker, Phys. Rev. B 66, 

024512 (2002); K. V. Samokhin and M. B. Walker, Phys. 

Rev. B 66, 174501 (2002); K. V. Samokhin, Phys. Rev. 

B 66, 212509 (2002). 



[4] I. A. Fomin, JETP Lett. 74, 111 (2001) [Pis'ma Zh. Eksp. 

Teor. Fiz. 74, 116 (2001)]; I. A. Fomin, JETP 95, 940 

(2002) [Zh. Eksp. Teor. Fiz. 122, 1089 (2002)]. 
[5] V. P. Mineev, Phys. Rev. B 66, 134504 (2002); V. P. 

Mineev, T. Champel, preprint cond- mat/0306471 to be 

published in Phys. Rev. B. 
[6] V. P. Mineev and K. V. Samokhin, Introduction to 



8 



Unconventional Superconductivity (Gordon and Breach, 
London, 1999). 

[7] P. W. Anderson and E. I. Blount, Phys. Rev. Lett. 14, 
217 (1965). 

[8] A. C. Lawson and W. H. Zachariasen, Phys. Lett. 38A, 
1 (1972). 

[9] J. F. Schooley, W. R. Hosier, and M. L. Cohen, Phys. 

Rev. Lett. 12, 474 (1964). 
[10] V. V. Bogatko and Yu. N. Venevtsev, Sov. Phys. Solid 

State 25, 859 (1983) [Fiz. Tverd. Tela 25, 1495 (1983)]. 
[11] V. M. Edelstein, JETP 68, 1244 (1989) [Zh. Eksp. Teor. 

Fiz. 95, 2151 (1989)]; V. M. Edelstein, Phys. Rev. Lett. 

75, 2004 (1995). 
[12] L. P. Gor'kov and E. I. Rashba, Phys. Rev. Lett. 87, 

037004 (2001). 

[13] V. Barzykin and L. P. Gor'kov, Phys. Rev. Lett. 89, 
227002 (2002). 

[14] L. S. Levitov, Yu. V. Nazarov, and G. M. Ehashberg, 

JETP Lett. 41, 445 (1985) [Pis'ma Zh. Eksp. Teor. Fiz. 

41, 365 (1985)]. 
[15] V. P. Mineev and K. V. Samokhin, JETP 78, 401 (1994) 

[Zh. Eksp. Teor. Fiz. 105, 747 (1994)]. 
[16] E. Bauer, G. Hilscher, H. Michor, Ch. Paul, E. W. 

Scheldt, A. Gribanov, Yu. Seropegin, H. Noel, M. Sigrist, 

and P. Rogl, Phys. Rev. Lett. 92, 027003 (2004). 
[17] V. P. Mineev, JETP Lett. 57, 680 (1993) [Pis'ma Zh. 

Eksp. Teor. Fiz. 57, 659 (1993)]. 
[18] E. Wigner, Group Theory and Its Applications to the 

Quantum Mechanics of Atomic Spectra (Academic Press, 

NY, 1959); C. J. Bradley and B. L. Davies, Rev. Mod. 

Phys. 40, 359 (1968). 
[19] L. D. Landau and E. M. Lifshitz, Quantum Mechanics, 

3rd Edition (Butterworth-Heinemann, Oxford, 2002). 
[20] E. I. Rashba, Sov. Phys. Solid State 2, 1109 (1960). 
[21] P. Blaha, K. Schwarz, P. Sorantin, and S. B. Trickey, 

Comp. Phys. Commun. 59, 399 (1990). 
[22] G. K. H. Madsen, P. Blaha, K. Schwarz, E. Sjostedt, and 

L. Nordstrom, Phys. Rev. B 64, 195134 (2001). 
[23] K. Schwarz, P. Blaha, and G. K. H. Madsen, Comp. Phys. 

Commun. 147, 71 (2002). 
[24] W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 
[25] J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 

(1992). 

[26] Some trial calculations using the linear muffin-tin orbital 
method [O. K. Andersen, Phys. Rev. B 12, 3060 (1975)] 
produced similar results for the LDA and the generalized 
gradient approximation [J. P. Perdew, J. A. Chevary, S. 
H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh, 
and C. Fiolhais, Phys. Rev. B 46, 6671 (1992)]. 

[27] A. H. MacDonald, W. E. Pickett, D. D. Koelling, J. Phys. 
C: Solid St. Phys. 13, 2675 (1980). 

[28] B. S. Chandrasekhar, Appl. Phys. Lett. 1, 7 (1962); 
A. M. Clogston, Phys. Rev. Lett. 9, 266 (1962). 

[29] A. I. Larkin and Yu. N. Ovchinnikov, Sov. Phys. JETP 



20, 762 (1965) [Zh. Eksp. Teor. Fiz. 47, 1136 (1964)]; 
P. Fulde and R. A. Ferrell, Phys. Rev. 135, 550 (1964). 
[30] That there is only one type of pairing possible in non- 
centrosymmetric systems with spin-orbit coupling, was 
pointed out in P. W. Anderson, Phys. Rev. B 30, 4000 
(1984). 

[31] P. Frigeri, D. F. Agterberg, A. Koga, and M. Sigrist, 
preprint |cond-mat /d3il354| (unpublished) . 



0.5 




-2 I ■ I ■ I ■ [• ■ 1 — I ' - I 

rXM rZRAM 

k vector 



FIG. 1; (Color online) Band structure calculated without SO 
coupling. Three bands (labeled a, /3, and 7) cross the Fermi 
energy. 
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FIG. 2: (Color online) Band structure of CePtaSi with SO 
coupling. Bands are split due to the SO interaction (cf. Fig. 
0. The bands that cross the Fermi energy, are plotted as 
solid lines. 
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FIG. 3: (Color online) Fermi surface of CePtsSi. The sub- 
figures (a)-(g) show cross sections for different values of kz 
{kz is measured in units of 2tt/c, and kx,y are measured in 
units of 2-n/a). Only one quarter of the Brillouin zone is 
shown. By applying the reflection symmetries in the fcy = 
and fcz = planes, the Fermi surface in the entire Brillouin 
zone can be recovered. The labels of the sheets correspond to 
the labels in Fig. Dots indicate the k points for which we 
calculated the band energies. The Fermi surface was obtained 
by interpolating between these points. 



